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ABSTRACT 

Context. Object with a very peculiar light-curve was discovered recently using Kepler data from first two quarters. Authors argue that 
this object may be a transiting disintegrating planet with a comet like dusty tail. Light-curves of some interacting binaries may have 
features analogous to this light-curve and it is very attractive to see if they are indeed caused by the same effects. 
Aims. The aim of the present paper is to verify the model suggested by the discoverers by the light-curve modelling and put constraints 
on the geometry of the dust region and various dust properties. 

Methods. We modify the code Shellspec designed for modelling of the interacting binaries to calculate the light-curves of stars with 
planets with comet like tails. We take into account the Mie absorption and scattering on spherical dust grains of various sizes assuming 
realistic dust opacities and phase functions and finite radius of the source of light (star). 

Results. The light-curve is reanalysed using first six quarters of the Kepler data. We prove that the peculiar light-curve of this objects 
is in agreement with the idea of a planet with a comet like tail. Light-curve has a prominent pre-transit brightening and a less prominent 
brightening after the transit. Both are caused by the forward scattering. Forward scattering is a strong function of the particle size. 
Dust density in the tail is a steep decreasing function of angle (distance) from the planet which indicates significant dust destruction 
along the tail caused by the star. The particle size of the grains in the tail is about 0. 1 micron but can be slightly larger if data with 
the shorter exposure (short cadence) were available. The orbital period of the planet was slightly improved. This light-curve with 
pre-transit brightening is anaologous to the light-curve of e Aur with mid-eclipse brightening and forward scattering plays significant 
role in both of them. 

Key words. Keywords should be given 



1. Introduction 

The exoplanet KIC0 12557548 has been discovered recently by 
Rappaport et al. (120121 1. It was discovered from Kepler Ql, Q2 
long cadence data (Borucki et al. 2011). This exoplanet is very 
unique. Unlike all other exoplanets it exhibits strong variability 
in the transit depth. For some period of time the transit even dis- 
appears. Shape of the transit is highly asymmetric with a sig- 
nificant brightening just before the eclipse, sharp ingress fol- 
lowed by a smooth egress. Planet has also extremely short pe- 
riod of 15.6854 hours. Rappaport et al. (2012 ) suggested that the 
planet has size not larger than Mercury and is slowly disintegrat- 
ing/evaporating what creates a comet like tail made of pyroxene 
grains. 

In this paper we revisit the light-curve using Kepler obser- 
vations from the first six quadratures and model the light-curve 
using the code Shellspec (Budaj & Richards 2004) which en- 
ables us to prove that the shape of the light-curve is in very good 
agreement with the idea of a planet with the comet like tail and 
put constraints on the particle size and geometry of the dust re- 
gion. 

2. Observations 

In this paper we use the publicly available Kepler data from the 
first six quarters (Borucki et al. 1201 II ). First we verify and im- 
prove the orbital period of the planet and search for periods in 
these data by two independent methods. Phase dispersion min- 
imization method (PDM, Stellingwerf 1978) was used first. We 



found two significant periods of 21.6 and 23.1 days which might 
be associated with the stellar differential rotation and some spots 
at different latitudes. However, our goal is to verify and improve 
the planet orbital period. Our first rough scan over that data indi- 
cates period of 0.653584 days which is very close to the period 
announced in the discovery paper. Consequently, we phased the 
data with this preliminary period and set the zero point such that 
the transits occur approximately at the phase of 0.5. Then we 
split the data into many chunks each covering one orbital period, 
fit the straight line into each chunk covering phases 0- 1 and re- 
move (divide out) the linear trend from the data. In this way we 
removed from the light-curve any changes on time scales longer 
than the orbital period. Then we again searched for the periods 
and found the final value of the frequency 1.530055 which cor- 
responds to 0.653571 days. This value is also in the good agree- 
ment with the value reported in the discovery paper. The error 
associated with the half-width of the spectral line correspond- 
ing to the lowest value of the theta is 0.0001 but our numeri- 
cal experiments with various parameters indicate that this error 
may be even 1 order of magnitude smaller. In the next step we 
used the Fourier method and found the frequency (period) of 
1.530063 eye (0.653568 days) which confirms the period found 
by the PDM method. 



3. Optical properties of the dust 



According to Rappaport et al. ( 120 121 ) a comet like tail of the 
planet mainly consists of dust maid of pyroxene grains. Optical 
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Fig. 1. Phase functions at 650 nm for different particle radii. Full 
line - phase functions assuming a point source of light. Dashed 
line -phase functions which take into account the finite dimen- 
sion of the stellar disk. Notice the strong forward scattering peak 
for larger particles. Vertical line illustrates the angular dimension 
of the stellar disk as seen from the planet. 
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Fig. 2. Light-curve with the transit. It is best represented by an 
exoplanet with a comet like tail. This tail is modeled as a part of 
ring - arc. Dust density in the tail rapidly decreases with the (an- 
gular) distance from the planet which indicates that they are ei- 
ther quickly destroyed or scattered out of the ring. Transit depth 
can be described with particles of different size and density. 



properties of the pyroxene: opacities for scattering and absorp- 
tion and phase functions were calculated using Mie theory and 
Mie scattering code of Kocifaj (120041 1 and Kocifaj et al. (|2008 ). 
The complex index of refraction of pyroxene was taken from 
Dorschner et al. (|19951 l. 

We assumed spherical particles of different sizes. We as- 
sumed a nonzero iron content where 20% of Mg atoms were 
replaced by Fe. To suppress the ripple structure that would ap- 
pear in the phase function of spherical mono-disperse particles, 
the poly-disperse Deirmendjian 119641 ) distribution of particle 
sizes was assumed. The phase functions for three different parti- 
cle sizes at 6500 A are illustrated in the Figure Q] 

However, these phase functions do not take into account the 
finite and non-negligible angular dimension of the stellar disk as 
seen from the planet. To take this effect into account one would 
have to split the stellar disk into elementary surfaces and inte- 
grate over the disk. Note that if the angular planet-star separa- 
tion is greater then the angular dimension of the star this process 
is mathematically analogous to the rotational broadening of the 
stellar spectra. Consequently, we convolve the phase functions 
with the 'rotational' broadening function (Rucinski 2002) of the 
star expressed in the units of the angular dimension of the star as 
seen from the planet. This function is an ellipse and takes into ac- 
count the linear limb darkening law. We used the limb darkening 
coefficients from Claret (2000) (see the next section). One has 
to be cautious and calculate the phase functions with a very fine 
step near the zero angle because of the strong forward scatter- 
ing and, consequently, the convolution with a very fine step near 
the edge of the broadening ellipse. These phase functions which 
take into account the finite dimension of the source of light are 
also depicted in the Figure Q] for comparison. Broadened phase 
functions are renormalized to 4-n during the calculations. 

4. Light-curve modelling 

In this section we calculate light-curves of the exoplanet us- 
ing the code Shellspec. The star is assumed to be a sphere 
with the radius of 0.65 R e , mass of 0.7 M and be subject to 
the limb darkening. We used quadratic limb darkening coeffi- 
cients for R filter from Claret (2000) which are based on the 



Kurucz models. They were interpolated from the grid assuming 
v, = 2km s~\ r eff = 4300^, logg = 4.63, [M/H] = 0. Given the 
fact that the transit is sometimes missing in the data the planet 
itself is rather small and it is its comet like tail which gives rise 
to the major drop in the light-curve. Consequently, the planet is 
modelled as a part of the ring - arc with the radius of 2.8 R Q . 
The cross-section, C, of the arc and dust density, p, along the arc 
is allowed to change with the angle, t in [rad], in the following 
way: 

Pit) = P (0)C(0)/c[|f - r(0)| + ir denrg . (i) 

Where p(0), C(0), f(0) is the dust density, cross-section and an- 
gle at the beginning of the arc, and edenrg is density exponent 
to model additional phenomena e.g. dust destruction. The best 
match (not necessarily a unique solution)^ with the observations 
is depicted in the Figure[2] It was obtained for the following tail 
parameters. The cross-section of the tail was a square which var- 
ied from 0.2x0.2 R at the beginning to 0.6x0.6 R e at the end. 
The light curve was not very sensitive to this cross-section. Tail 
was 60 degrees long and edenrg m -7. It means that it expe- 
riences a considerable dust destruction along the tail. As men- 
tioned in the discovery paper this weakening tail is responsible 
for the extended absorption in the light-curve. 

A close look at the light-curve reveals a significant pre- 
transit brightening and a weak brightening after the transit. It is 
very interesting to note that this pre-transit brightening is anal- 
ogous to the mid-eclipse brightening observed in some eclips- 
ing interacting binaries, particularly e Aur (Budaj 201 lb). Mid- 
eclipse brightening in e Aur may be caused by the flared dusty 
disk and forward scattering plays an important role. So although 
these are very different objects they are subject to similar physi- 
cal processes and can be modelled with similar tools. As demon- 
strated in FigureQ~l forward scattering is strongly sensitive to the 
particle size and, consequently, we used this feature to estimate 
that dust particles have radii of about 0. 1 micron. See Figure[3] 



1 I tried also a number of different models (tail of the shape of a cone, 
exponential density profile along the tail) but the present model seems 
to describe the observations better. 
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Fig. 3. A detail look at the light-curve. Pre-transit brightening is 
caused by the strong forward scattering. It is rather sensitive to 
the particle size which can be constrained reasonably well. Dust 
particles in the tail have typical radii of about 0.1 micron. 



5. Conclusions 

The main results are the following: The light-curve of this planet 
was reanalysed using first six quarters of the Kepler data. Orbital 
period of the planet was improved. We modelled the light-curve 
of KIC0 12557548 using Shellspec code assuming realistic dust 
properties (pyroxene spherical grains of different particle sizes, 
Mie absorption and scattering) finite radius of the source of light 
(star) and that the medium is optically thin. We proved that its 
peculiar light-curve is in agreement with the idea of a planet 
with a comet like tail. Light-curve has a prominent pre-transits 
brightening and less prominent brightening after the transit both 
caused by the forward scattering. Dust density in the tail is a 
steep decreasing function of angle (distance) from the planet 
which indicates significant dust destruction along the tail caused 
by the star. Transit depth is a highly degenerate function of the 
particle size and other dust properties. However, forward scatter- 
ing is strongly sensitive to the particle size. We estimated the par- 
ticle size of the grains in the tail from the pre-transit brightening 
(0.1 micron). It is interesting to note that this planet light-curve 
with pre-transit brightening is analogous to light-curves of some 
interacting binraries with mid-eclipse brightening, particularly e 
Aur (Budaj 2011b) and forward scattering plays important role 
in both of them. 
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2 I realized that very recently similar calculations and conclusion 
were obtained independently by Brogi et al. d2012l >. These authors used 
analytical Henyey-Geenstein phase functions. 
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